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We calculate the phase diagram for dilute mixtures of bosons and fermions at zero temperature. 
The linear stability conditions are derived and related to the effective boson-induced interaction 
between the fermions. We show that in equilibrium there are three possibilities: a) a single uniform 
phase, b) a purely fermionic phase coexisting with a purely bosonic one and c) a purely fermionic 
phase coexisting with a mixed phase. 



I. INTRODUCTION 

Recent developments in the trapping and cooling of 
atoms have made it possible to investigate the proper- 
ties of dilute gases at very low temperatures where the 
bosonic or fermionic character of the atoms becomes cru- 
cial. Following the realization of Bose-Einstein condensa- 
tion in a single-component gas of alkali atoms fl]-|| , the 
group at JILA succeeded in trapping and cooling 87 Rb 
atoms in two different hyperfine states, thereby creating 
overlapping condensates in this boson-boson mixture [Q . 

Lately a lot of attention has been given also to 
fermions, especially in view of the possibility of achieving 
temperatures low enough to observe a BCS transition. 
Due to the vanishing s-wave cross section for identical 
fermions in the same spin state, evaporative cooling of 
a single species of fermions is ineffective. By using a 
mixture of two spin states DeMarco and Jin || recently 
succeeded in cooling the fermionic isotope 40 K to temper- 
atures lower than the degeneracy temperature. Another 
possibility, however, is to trap fermions together with 
bosons and cool the latter, so that the fermions are cooled 
through their thermal contact with the bosons (so-called 
sympathetic cooling). Such fermion-boson mixtures form 
the topic of the present paper. 

The miscibility properties of boson-boson mixtures in a 
trap have been discussed in several papers § §. Trapped 
boson-fermion mixtures have been considered in ref . |9f | 
within the Thomas-Fermi approximation and in ref. [ JLOf , 
where the separation of the components was studied nu- 
merically as a function of the interparticle interaction. 
The purpose of the present work is to carry out an ana- 
lytical study of the miscibility of fermion-boson mixtures 
in the uniform case, in order to obtain insight into the 
types of phase boundaries that may occur in this system. 

The article is organized as follows. We begin in sec- 
tion [n] by studying linear stability of uniform mixtures 
and introduce the effective fermion-fermion interaction 
as mediated by the bosons. In section III wc determine 



the general conditions for phase equilibrium by consider- 
ing the pressure and chemical potentials on each side of 
a phase boundary. We find that depending on the total 
densities of the two components, as shown in Fig. 3, there 



are three possibilities, a uniform phase with bosons and 
fermions fully mixed, a purely fermionic phase coexisting 
with a purely bosonic one, and a purely fermionic phase 
coexisting with a mixed phase. However, a separation 
into two phases, each with different, non-zero, concen- 
trations of bosons and fermions is never in equilibrium, 
nor is one with a purely bosonic phase and a mixed one. 
In section IV we discuss the experimental implications of 
our results, using the fermion isotopes 6 Li, 40 K and 84 Rb 
as examples. 



II. LINEAR STABILITY AND INDUCED 
INTERACTIONS 

The general conditions for stability of a binary mixture 
towards small changes in the concentrations of its com- 
ponents are obtained by considering the total energy, £ , 
as a functional of the densities n a and n$ of the two 
components: 



£ = drE{n a (r),nf3(r)). 



(1) 



We consider the change in total energy arising from small 
changes, 5n a and Snp, in the concentrations of the two 
constituents. The first order variation, S£ , must vanish, 
since the number of particles of each species is conserved, 



drSrii = 0; i = a,j3. 



(2) 



The second order variation, 5 2 £ , is given by the quadratic 
form 
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In carrying out this expansion we have implicitly assumed 
that the characteristic wavelength for the spatial varia- 
tion of the densities is long compared to the microscopic 
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lengths in the problem. As we shall see, for dilute mix- 
tures of bosons and fermions the relevant length is the 
coherence length of the bosons. The derivative of the en- 
ergy density with respect to the particle density, dE/drn, 
is the chemical potential for the species labelled i 
(i = a, 0). The quadratic form (||) is thus positive semi- 
definite, provided 



dn a 



>0, 



9fJ£ 

dn g 



>0, 



and 



dn a drip dnp dn a 



> 0. 



(4) 



(5) 



Note that it is only necessary that one of the two condi- 
tions in (||) is satisfied, since the other one is automati- 
cally fulfilled when the condition (|5|) is. 

For a dilute mixture of bosons with fermions in a single 
internal state, at zero temperature the energy functional 
has the form 



£ = V 



n B n F U B F + -= e FnF ) 



(6) 



where V is the total volume. The first two terms are 
due to the boson-boson and boson-fermion interactions, 
respectively, with n B and np denoting the boson and 
fermion densities. Ubb is the matrix element of the ef- 
fective interaction for bosons with bosons and Ubf that 
for bosons with fermions. The last term represents the 
kinetic energy of the fermions. The fermion- fermion in- 
teraction energy is negligible, since a) the s-wave scat- 
tering amplitude vanishes for fermions in the same spin 
state and b) none of the higher partial waves contribute 
at the low temperatures we are considering. The chemi- 
cal potentials are then seen to be 



fiB = n B UsB + tifUbf, 
Hf = £f + u b Ubf- 



(7) 



Let us introduce the constant A through the definition 



where tub is the boson mass and tubf — 2m f^b/ (m f + 
mg) is twice the reduced mass of a boson-fermion pair, 
then the condition (||) becomes 
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If the masses as well as the scattering lengths occurring 
in ( pd| ) are approximately equal in magnitude, then the 
stability condition requires the mean interfermion dis- 
tance, which is approximately rip 1 ^ , to be greater than 
a scattering length. The trapped one-component gases 
that have so far been investigated experimentally are di- 
lute in the sense that the mean interparticle distance is 
much greater than the scattering length. The stability 
condition ( pi] ) would therefore generally be expected to 
hold for dilute mixtures, unless the scattering length clbf 
greatly exceeds cibb- In this case the condition ( pd| ) 
may be violated at much lower fermion densities given 



by n)/ 3 aBB ~ (o-bb/o.bf) 2 rather than rip^asB ~ 1- 

It should be noted that the mean interfermion distance 
in the stability condition ( pi] ) formally plays the role of 
an effective fermion-fermion scattering length. This is 
apparent when we consider binary mixtures of bosons, 
labelled by a and (3. In this case the energy functional 
takes the form 

£ = V (~n 2 a U aa + n a n p U a (3 + ^npUpp^] • (12) 



The condition (^) would then be replaced by 



(13) 



The conditions (^) and (|l3|) assume the same form, 
apart from numerical factors, provided the combination 

AirT^rip 1 ^ /mp is interpreted as arising from a direct 
fermion-fermion interaction, while physically, of course, 
the term originates in the kinetic energy. The Fermi pres- 
sure thus behaves as a two-body interaction with a scat- 
tering length of order the particle separation. 
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where ttif denotes the fermion mass. The value of the 
numerical factor in A reflects the fact that we consider 
fermions in a definite spin state. According to (|^) linear 
stability then requires 
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U BF 



(9) 



If the interaction parameters are converted to scattering 
lengths asB and asF through the relations 

4irh 2 a B B TT Anh 2 a B F MnS 
Ubb = ; Ubf = , (1U) 



m B 
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A. Induced interactions 

In a multi-component system the effective interaction 
between atoms of one species is altered by the pres- 
ence of the other components. As an example of this 
we return to the stability condition (|J). Since the bo- 
son gas must be stable to density fluctuations when the 
fermions are uniform, the compressibility of the bosons 
must be positive, which means that d[isjdnB > 0. Also, 
diiB/dnp = d[iF/dn,B = d 2 E/dnBdnp. Thus the sta- 
bility condition (g|) may be written as 
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(14) 
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The first term here is the variation of the chemical poten- 
tial of the fermions when the fermion density is changed, 
keeping the boson density fixed. The second term may be 
regarded as an induced interaction, due to the fact that 
one fermion tends to attract or repel bosons, depend- 
ing on the sign of the boson-fermion interaction, and the 
change in the boson density changes the energy of a sec- 
ond fermion. The induced interaction is analogous to the 
phonon-induced attraction in metals. Note that this in- 
teraction is always attractive, irrespective of the sign of 
the boson-fermion interaction. By use of the chain rule, 
one may show that the stability condition is equivalent 
to 



driF 



> 0. 



(15) 



This way of looking at interactions has been ex ploi ted in 
the context of liquid mixtures of 3 He and 4 He |X1| . The 
induced interaction contribution to the stability criterion 
corresponds to the diagram shown in Fig. 1, evaluated in 
the long- wavelength limit. 



U ind {q,u}) = U BF XB{q,u). 



(19) 



In the Bogoliubov approximation the density-density re- 
sponse function is given by 



n B q 



(20) 
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where the excitation frequencies are the Bogoliubov ones 
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(21) 



with e q = h q /2m B being the free boson energy. Thus 
the static induced interaction is 

U in d(q,oj) = -U BF 



n B U BB + (h' ! q 2 /±m B ) 



(22) 



In coordinate space this is a Yukawa, or screened 
Coulomb, interaction 



^ind(r) = - 



m B n B U BF e 



-V2r/£ 



(23) 




F / \ F 

FIG. 1. Diagrammatic representation of the interaction be- 
tween fermions induced by exchange of a boson density fluc- 
tuation. 



For a dilute gas, the induced interaction is easily eval- 
uated using the expression (§) for the energy, since 



9hf 
dn B 



= U BF 



and 



d[i B 
dn B 



= U BB . 



Thus the induced interaction is given by 

J I 2 

TT — BF 

Ubb 



(16) 



(17) 



(18) 



and it is equivalent to an effective reduction of the scat- 
tering length by an amount {a 2 BF / a BB ){m B mp / m 2 BF ) . 

At non-zero wavenumbers and frequencies the effective 
interaction is the natural generalization of Eq. ( |l8| ) , the 
static long- wavelength response function —dn B /d/j, B for 
the bosons being replaced by the density-density response 
function \b (q, w). Thus 



where £ is the coherence (healing) length for the bosons, 
given by 



2m B n B U BB 



(24) 



The induced interaction thus has the interesting fea- 
ture that at long wavelengths it gives rise to effects which 
are of the same order of magnitude as a typical bare in- 
teraction, if the boson-boson and boson-fermion interac- 
tions are of the same order of magnitude and boson and 
fermion masses are comparable. The reason for this is 
that even though the induced interaction involves two 
boson-fermion interactions, the density-density response 
function for the bosons at long wavelengths is inversely 
proportional to the boson-boson interaction. At wave 
numbers greater than the inverse of the coherence length 
for the bosons, the magnitude of the induced interac- 
tion is reduced, since the boson density-density response 
function has a magnitude ~ 2n B /e q . The induced inter- 
action is thus strongest for momentum transfers less than 
m B s B , where sb — (tibUbb /wb) 1 ' 2 is the sound speed 
in the boson gas. If bosons and fermions have comparable 
masses and densities, and the scattering lengths are com- 
parable, for momentum transfers of order the Fermi mo- 
mentum the interaction will be of order A-Kh 2 k F a 2 j '(2m). 
The induced interaction is therefore one power of the 
"diluteness parameter" k F a less than a typical direct in- 
teraction. For mixtures of two species of fermions (for 
example two different hyperfine states) with bosons the 
induced interaction is attractive, and therefore it will in- 
crease the transition temperature to a BCS superfluid 
state. Detailed calculations are necessary to determine 
how important the effect is quantitatively. 
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III. PHASE EQUILIBRIUM 

The stability considerations given in the previous sec- 
tion are valid only for small changes. We now solve the 
general stability problem by considering the conditions 
under which mixtures with different concentrations of 
bosons and fermions may be in equilibrium with each 
other. 

We consider a mixture of Ng bosons and Np fermions 
in a box of volume V. These components may either mix 
uniformly or form distinct phases, which we label by the 
index i. If we ignore interpenetration effects, a possible 
phase-separated configuration may be described by the 
number of phases present, /, the bosonic and fermionic 
densities in each phase, tiba and np.i, and the fractions of 
the total volume they occupy, Uj. Since the total number 
of particles is given, the following relations must hold: 
X«=i n F,iV t = n F = N F /V, Yli=i n BAVi =n B = N B /V, 
and v i = 1- When I = 1 we recover the case of a 
homogeneous mixture. Let us now turn our attention to 
the case 7 = 2. The total energy is the sum of contri- 
butions due to boson-boson interactions, boson- fermion 
interactions and the kinetic energy of the fermions, 

2 2 

£ = J2£ i = Vj2v i E i , (25) 

»=1 i=l 

where 

1 3 

Ei = 2 n2 B,i U BB + n B ,i n F,iU BF + ^ F ^n F<i , (26) 

with the Fermi energies e F j being given by 

%2 

la 2n2/3 2/3 n / 07 x 

e Fji = (6n ) n^.—. (27) 

The pressure pi in each phase is then found by differ- 
entiating the energy with respect to the volume V$ = Vvi 
occupied by each phase, 

dSi 1 2 2 

Pi = = 2 U B,i U BB + n B ,i n F,iU BF + ^ F ^n F>i , 

(28) 

usingthe volume dependence of the Fermi energy given 
by (|27|). The first requirement for equilibrium is that 
Pi = P2- The two other conditions involve the chemical 
potentials which are given by 

dE % 

(J-B.t = q = nB,iUBB + n F iUBF, 

oriB,i 
dEi 

Ufa = q — - = CF.j + n B ,iU B F- (29) 
on F ,i 

If the boson density Uba is non-zero in both phases, then 
the chemical potentials (jlba must be equal. If the boson 
density vanishes in one phase, then the boson chemical 



potential in that phase must be higher than in the other 
one, in order for the system to be in equilibrium. The 
same is of course true for the fermions. 

We remark that for a two-component system it is not 
possible to have 3 or more distinct phases in equilibrium. 
This is seen as follows. Consider a situation where there 
are I distinct mixed phases. The total number of condi- 
tions to be met is 3(7 — 1). If one wishes equilibrium to 
be possible over a range of parameters one may in addi- 
tion fix another variable, such as the total pressure. The 
total number of conditions to be satisfied is then 31 — 2. 
However the number of independent variables is just 21, 
corresponding to the boson and fermion densities for the 
I phases. It is then clear that solutions to the equations 
are only possible if I < 2. Similar arguments apply if 
some pure phases are present. 

We thus consider a system consisting of two phases 
with densities tib,i,ti Fj i and riB,2,nF,2, with the first one 
occupying a share v of the total volume. Due to the con- 
ditions on the chemical potentials one has to distinguish 
four cases, which must be analyzed one by one: 

A. Two pure phases: The bosons and fermions are com- 
pletely separated corresponding to n F i = 0, rig, 2 = 
and riB.i ^ 0, n F ^ 0. 

B. A mixed phase and a purely fermionic one: The boson 
density vanishes in one region corresponding to n^i ^ 0, 
ns,2 = and n B s ^ 0, n F2 ^ 0. 

C. A mixed phase and a purely bosonic one: The fermion 
density vanishes in one region corresponding to the con- 
ditions n F _i = 0, Ub,2 7^ and fi^.i 7^ 0, n Fj 2 =/= 0. 

D. Two mixed phases: All densities tib,i, np^i and ns,2, 
rip, 2 arc different from zero, while ns.i ns.2 and upj ^ 
n Ft 2- 

First we consider cases A and B and identify the two 
regions of the ns-np plane in which they can occur. We 
demonstrate below that these regions are in fact different 
and that in each of them the phase-separated configura- 
tion is energetically favored compared to the homoge- 
neous one. Subsequently we prove that the cases C and 
D are not realizable. 



A. Two pure phases 

As mentioned previously the pressures in the two 
phases must be the same. These are found from the gen- 
eral formula (|28| ) by setting n F ^\ — in the expression for 
Pi and riB,2 = in the expression for pi. The condition 
Pi = P2 then becomes 

^n 2 B1 U B B = -An^l, (30) 

which shows that there is only one allowed value for ub,i 
for a given fermion density n F ,2- 

Let us now look at the chemical potentials. If the 
chemical potential in the fermionic phase 2 has to be 
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lower or at most equal to that in the bosonic phase 1, 
then it follows from ( p9|) that 



£F,2 < riB.lUBB- 



Likewise, if the chemical potential in the bosonic phase 
has to be lower than in the fcrmionic one, then 



riB,iUBB < np^UsF, 



(32) 



which shows that Ubf must be positive. Inserting the 
value of ns.i from (^0|) into and ( |32"|) we obtain from 
( |3l| ) the condition 



n F ,2 > 
while @ yields 

n F ,2 > 



5AU. 



( 4f? 



BB 



BF 



( 4AU 



BB 



BF 



(33) 



(34) 



which is less restrictive than (|33|). Hence equilibrium of 
a pure boson phase and a pure fermion one is possible if 
and only if the fermion density satisfies (|33| ) and ub,i is 
given by @. 
Now, since 



tlb = v riB,i and 
n F = (1 - v) n F ,2, 



(35) 



where v is the fraction of the total volume occupied by 
the bosons, we can find those total densities for which 
the system can completely separate by using the allowed 
values of n F ,2 and Ub,\ in (|35|). These densities corre- 
spond to the points above the uppermost full line in the 
rig vs. n F plot in Fig. 3. 

Let us now compare the energy of a phase-separated 
configuration with that of a mixed state with the same 
number of particles, as given by (|3^). According to (||) 
the energy of the mixed system is 



f 



V 



1 3 N 

-n B U BB + n B n F U BF + -n F e F 



(36) 



where the densities n B and n F are given by (|35|). On the 
other hand, the energy of the phase-separated system 
with the same number of particles is 



1 3 
£ SC p = V [ -n B1 vU B B + -eF,2"F,2(l 



(37) 



The energy difference £ m ix — £scp is obtained by substi- 
tuting for us and n F from (Bq), with n B \ given by (30) 



-(£mix - £scp) = - V ) An F,l 



V 



+ v(l - v^^Ubp^A/Wbb) 112 

- |(1 -v)(l -(l-vf/ 3 )An 5 Jl (38) 



Since by ( ^3| ) we must have n)/® > \/5AUbb/4U BFi 
then 



( 31 ) i(S mix - £ scp ) > 3 -An% 3 2 (l - v) 5 ^(l (1 - «) 1 /3) > o 



(39) 



for any v between and 1, and the equality holds only at 
V = and v = 1 . So in the density ranges where equilib- 
rium is possible £ scp is always less than £ m ; x , and phase 
separation is energetically favored compared to mixing. 



B. A mixed phase and a purely fermionic one 

We now consider a mixed phase in equilibrium with a 
pure fermion phase and let ub,2 = 0. As we shall see in 
the following, this type of phase separation is stable in 
a different range of total densities than that for case A. 
The pressure equilibrium condition takes the form 



1 



\Ubb + n B ,in Ftl UBF + -£ F ,in F ,i = -e F ^n Ft 2- 



(40) 

Beside this we have to impose also the equality of the 
fermion chemical potentials in the two phases, which im- 
plies 



A 
Ub^ 



I 2/3 
\ n F2 



2/3 



(41) 



First let us assume that Ubf is positive. Then ( |4l| ) 
implies that uf,2 > ^f,i, and it allows us to eliminate 
Ub,i from equation (40). The resulting equation is con- 
veniently expressed in terms of the ratio between the 
fermion densities, 



»F,1 

n F ,2 



1/3 



and the dimensionless constant A, defined by 



A = 



AU 



BB 



U BF ' 



l 1/3 ' 
l F ,2 



Inserting hb.i from ( pd] ) into ( [40| ) we thus obtain 



A 



(42) 



(43) 



(44) 



By solving this equation one obtains as a function of 
n,f t 2- The problem is greatly simplified by noticing that 
(1 — x) 2 is a common factor on both sides of Eq. (Q). The 
value x = 1 implies that the fermion density is the same 
everywhere, while the boson density vanishes everywhere 
according to eq. (J4l| ). The solution x — 1 can therefore 
be rejected, and ( f44| ) becomes the cubic equation 

6x 3 + (12 - 5A)a; 2 + (8 - I0\)x + 4 - 5A = (45) 
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Since we have assumed that Ubf is positive or zero, 
which, according to (fll|), implies that np.i < ^f,2> the 
physical range of x lies between and 1. Evidently 
A = 4/5 yields the solution a; = 0, while i = 1 is a 
solution provided A = 3/2. The cubic equation ( f45| ) has 
one positive root in the range < x < 1, when A lies in 
the range 



4 3 
5 <A< 2- 



(46) 



According to ( [43|) this implies that the fermion density 
ni?2 is bounded from below as well as from above. We 
conclude that to each value of uf,2 lying within the range 
(46), there corresponds one value of np,i and, from Eq. 
(41), one value of ub.i- 

In Fig. 2 we plot the values of np,i and Ub,i found for 
each HF2- 



« 0.5 




U F,2 

FIG. 2. Densities of fermions (tif,i) and bosons (ns,i) in a 
mixed phase in equilibrium with a pure fermion phase of den- 
sity 71f,2. The fermion densities are in units of A 3 U bb /U B f, 
and the boson one in units of A i U% B /U% F . The vertical line 
corresponds to the value (2/3) 3 , for which uf,i = «f,2 and 
ub.i = 0. The dashed parts of the curves are unphysical for 
Ubf > 0, since the boson density is negative. For negative 
Ubf they correspond to solutions with positive boson den- 
sity, since the boson densities in the figure are then scaled to 
a negative quantity. However for the mixed phase under these 
conditions n^i > (2/3) i A i U BB /U%F an d it is unstable. 



us = v ns,i and 

n F = v uf,i + (1 — v) np,2, 



(47) 



they are obtained by using allowed values for np^, Uf,i 
and tib,i, and varying v from zero to one. These densi- 
ties correspond to the points in the region between the 
upper and lower full curves in Fig. 3. The lower line is 
obtained by letting v = 1 and is therefore just the set of 
allowed ns.i and uf,i- The upper one is found by set- 
ting np :2 = (5/4) 3 (AUbb/U bf ) 3 , the maximum value 
according to ([49), to which correspond uf,i — and 
ns.i = (5/4) 2 {A 3 U BB /U BF ), and varying v. It should 
be noted that this upper curve is exactly the one above 
which complete separation into a pure boson phase and a 
pure fermion one could take place. This proves the state- 
ment that cases A and B are stable for different values of 
total densities, which in view of the proof on the energy 
of case A, also demonstrates that above the upper line of 
Fig. 3 the system will be completely phase separated. 

Also in the present case we can show that the phase- 
separated configuration has lower energy than the cor- 
responding mixed one. The proof proceeds as in case A 
and we shall therefore omit it. 
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FIG. 3. Phase diagram of uniform boson-fermion mix- 
tures. The total fermion density rip is given in units of 
A 3 Ugg/U B p, and the total boson density ub in units of 
A 3 U B b/Ubf- The horizontal line corresponds to the value 
n F = (2/3) 3 A 3 U B b/U B f> above which the uniform system is 
unstable to small fluctuations. 



When a mixture is separated into two phases with one 
allowed set of tif,2, n F,i and Ub,i, it is straightforward 
to show that the boson chemical potential is always less 
in the mixed phase than in the purely fermionic one, but 
for brevity we omit the proof. 

Just as we did in the previous section, we can deduce 
the set of all total densities which could undergo a phase 
separation of this type. Since 



One more comment is due at this point. In the linear 
stability analysis we found that the homogeneous system 
would become unstable at densities violating (||) , but this 
value of the fermion density is just where the lower curve 
in Fig. 3 meets the vertical axis. 

Finally let us consider the situation when Ubf is neg- 
ative. In this case the fermion density in phase 1 must be 
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grea ter than that in phase 2, np,i > 1*^,2, according to 
([il|). We find again that the density ratio x — np^i/np^ 
must satisfy ((45|), but now x lies in the region between 
1 and infinity. By solving the cubic equation (^5|) it is 
readily shown (compare Fig. 2) that the values obtained 
for Uf,i are now greater than that which is required to 
ensure linear stability in phase 1. The phase-separated 
configuration is thus unstable for Ubf < 0. The at- 
tractive boson-fermion interaction in phase 1 leads to a 
collapse of the configuration with a mixed and a purely 
fermionic phase, presumably to a state where the bosons 
and fermions are clumped together and where the total 
energy is no longer given by (|6j). 

With this our knowledge of the phase diagram is com- 
plete, since the cases C and D cannot be realized, as 
proven below. 



C. Two mixed phases 

We shall now demonstrate that the cases C and D can- 
not be realized. First we consider two phases with differ- 
ent non-zero boson and fermion densities. 

The condition for equality of pressures is 

o( n l,l _ "^B.^BB + (nB,lTlF,l - n B,2 n F,2)U BF 
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0. 



and for equality of the chemical potentials 

(«B,1 - n B ^)U B B = (nF,2 - n F ,i)U B F, 
for bosons and 

(«B,1 - n B ,2)U B F = A(n 2 Jl - n 2 Jl). 



for fermions. Combining equation (|49|) with 
using the definitions and (E3J) we get 



1 



A(l 



(48) 

(49) 

(50) 
and 

(51) 



Substituting (^9|) into ((48|) and rearranging terms, we 
obtain furthermore 



(1 



= -X(l-x% 



(52) 



By writing (1—x ) = (l + x 3 )(l — x 3 ) and inserting l — x 3 



from fell) into (p2J) we obtain 




0. 



(53) 



The left hand side of ( |53| ) has a triple root at x = 1, and 
(p3|) can therefore be written as 



-{x 2 +3x + l)(x - l) 3 







(54) 



which clearly shows, since x is positive, that the only 
solution is the trivial one, x = 1. We have thus proven 
that two mixed phases cannot be in equilibrium with each 
other. 



D. A mixed phase and a purely bosonic one 

Finally we consider the case when one of the phases, 
say phase 2, is free of fermions and prove that equilibrium 
is impossible. We observe that (^8|) and ( f49| ) are valid 
also in this case, if we put Uf,i — 0. The same is true 
also for equation (^2j) , derived by combining them, which 
upon the use of x = becomes 



A 



(55) 



According to the definition ( ^3[ ) of A the system can 
therefore only be in equilibrium for one particular value 
of the fermion density. However, the fermion chemical 
potential in region 1 must be greater than that in region 
2, 



nB,iU B F > Aripl 



(56) 



Substituting for ub,i — nB,2, as derived from equation 
([ilil) with np,i = 0, we obtain 



Ubf ^ a 2/3 
Ubb 



(57) 



which according to the definition ( f43| ) implies A < 1. 
But the solution ( |55| ) does not satisfy this condition and 
must therefore be discarded. We conclude that case D 
also cannot be realized physically. 



IV. DISCUSSION AND CONCLUSIONS 

As we have emphasized, the properties of boson- 
fermion mixtures are of interest because such mixtures 
allow sympathetic cooling of fermions, and because the 
boson-induced fermion-fermion interaction, being attrac- 
tive, favors pairing of fermions, thus raising the BCS 
transition temperature. In this concluding section we 
discuss the implications of our results for actual experi- 
mental conditions. 

Our main findings are summarized in Fig. 3. This 
phase diagram shows that fermions and bosons under 
most circumstances will form a mixed phase, unless the 
interaction parameters assume quite large values. We 
shall now quantify this statement, using lithium, potas- 
sium and rubidium atoms as concrete examples. From 
the results obtained above, for instance the stability con- 
dition ( |33"| ) , it is apparent that the value of the dimension- 
less fermion density np(U BF / 'AUbb) 3 plays a key role 
in determining whether the fermions and bosons mix or 
phase-separate into either a purely fermionic phase to- 
gether with a mixed phase or a purely fermionic together 
with a purely bosonic phase (cf. Fig. 3). The correspond- 
ing dimensionless value of the boson density is according 
to © given by n B (U BF /A 3 U 2 BB ). 

Let us therefore look at the values of the actual densi- 
ties for which these dimensionless densities are equal to 
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one. A mixture of Li (boson) and 6 Li (fermion) has re- 
cently been produced |12| with 7 Li trapped in the hyper- 
fine state |1, —1} which, contrary to the |2, 2) state used 
in ||, has a positive scattering length sa 10 a.u., see [ p^| . 
For a 7 Li- 6 Li collision, on the other hand, the scattering 
length is « 38 a.u. H|. From this we deduce that np = 1 
corresponds to the density np = 9.3 • 10 23 m~ 3 , which 
is several orders of magnitude higher than the densities 
which can presently be realized experimentally. When 
expressed in terms of the boson-boson scattering length 
asB the density is seen to be np — 1.2 • 10~ 4 aJ 3 
eating the smallness of the gas parameter npa BB . 

Other interesting mixtures involve the potassium iso- 
topes 39 K, 40 K and 41 K, of which the first and the last are [4] 
bosons. In these cases using the values for the scattering 
lengths given in ]l5fl , we predict that phase separation ef- 
fects would take place for densities exceeding a B 3 a , which 



*a BB , indi- 



"BBi 

not only are much higher than can be experimentally re- 
alized, but also fall outside the range of applicability of 
the dilute gas approximation. 

A possible candidate for phase separation does, how- 
ever, exist among the rubidium isotopes. From the scat- 
tering lengths given in Jl6j we find that for 87 Rb and 
84 Rb mixtures the phase separation could take place at 
reasonable densities. Using the estimates asB ~ 100 
a.u. and cibf ~ 550 a.u., we have np = 1 when np = 
8.1 • 10 _6 a iJ 3 g = 7.1 ■ 10 19 m~ 3 , which is of the same order 
as the densities already realized in experiments with Rb 
isotopes. 

It is also illuminating to compare our findings for the 
uniform case with previous calculations |^,|l^| for mix- 
tures in a trap. Due to the Fermi pressure, the fermion 
density in a trap is usually much smaller than the bo- 
son density, and it is therefore a good approximation 
to assume that the bosonic cloud is unaffected by the 
fermions. The fermions on the other hand see an effec- 
tive potential, which may be attractive or repulsive inside 
the boson cloud, depending on the potential parameters 
and the ratio Ubf/Ubb of the interaction parameters 
If for simplicity the two trapping potentials are 
taken to be isotropic with same force constants, the ef- 
fective potential is repulsive if Ubf > Ubb- In this case 
the fermions are therefore expelled and form a shell out- 
side the boson condensate. If the effective potential is 
attractive, corresponding to the condition Ubf < Ubb, 
the fermions will reside inside the boson cloud. From 
the discussion above it is clear that in most cases the 
fermions inside the cloud will be fully mixed with the 
bosons. If the interaction parameters are appropriately 
chosen one then enters the phase-separated regime of our 
phase diagram, as is the case in the example shown in 
Fig. 2 of |9j, where a pure fermion phase forms at the 
center of the trap. The phase diagram which we have 
obtained here for the uniform case may thus serve as a 
useful guide for understanding the general behaviour of 
trapped mixtures. 
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